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Following the discovery of superconductivity in 
an iron-based arsenide LaOi j;Fj;FeAs with a su- 
perconducting transition temperature (Tc) of 26 
K[l], Tc was pushed up surprisingly to above 
40 K by either applying pressure [2] or replac- 
ing La with Sm[3], Ce[4], Nd[5] and Pr[6]. The 
maximum Tc has climbed to 55 K, observed in 
SmOi^F;j;FeAs[7|, HI and SmFeAsOi ;j;[9]. The 
value of Tc was found to increase with decreasing 
lattice parameters in LnFeAsOi j;Fj; (Ln stands 
for the lanthanide elements) at an apparently op- 
timal doping level. However, the F~ do ping in 
GdFeAsO is particularly difficult [l^, [ill due to 
the lattice mismatch between the Gd202 layers 
and Fe2As2 layers. Here we report observation of 
superconductivity with Tc as high as 56 K by the 
Th''+ substitution for Gd^+ in GdFeAsO. The in- 
corporation of relatively large Th^+ ions relaxes 
the lattice mismatch, hence induces the high tem- 
perature superconductivity. 

LnPeAsO family^] crystaUine in tetragonal 
ZrCuSiAs-typeflJl structure with space group P4/nmm. 
From the crystal chemistry point of view, the crystal 
structure can be described as an alternate stacking of 
Ln202 fluorite-type block layers and Fe2As2 antifluorite- 
type block layers along c-axis. The two block layers are 
connected by CsCl-type layers (Figure 1). Therefore, 
the chemical stability of LnFeAsO depends, to some 
extent, on the lattice match between the two block 
layers. A rough estimate based on the effective ionic 
radii [14] indicates that the Ln202 planar lattice is 
substantially smaller than the Fe2As2 lattice. The 
lattice mismatch becomes more serious for the LnFeAsO 
member with a smaller Ln^"*" ion. That may explain 
why TbFeAsO, DyFeAsO and other heavy- lanthanide- 
containing LnFeAsO members were hard to synthesize 
previously |12|. 

As the family member with the relatively small Ln'^+ 
ion, GdFeAsO is a promising parent compound to have a 
higher Tc through carrier doping. Similar to cuprate su- 
perconductors in which superconductivity emerges when 
charge carriers are induced into Gu02 planes by chem- 
ical doping at "charge reservoir layers", 4 5] supercon- 
ductivity in LnFeAsOi_a;Fa; is realized by partial sub- 
stitution of O^" with F-. The F^-for-O^" substitu- 
tion introduces extra positive charges in the insulat- 
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FIG. 1: Crystal chemistry understanding of the struc- 
ture of LnFeAsO (Ln=lanthanides). The stacking of flu- 
orite (CaF2) layers, CsCl-type layers and antifluorite (Li20) 
layers along c-axis forms the LnFeAsO structure. The lattice 
constant along the stacking direction can be expressed by the 
formula c ~ ^acaF2 + \acsCi + ffflLiaO, which basically sat- 
isfies the experimental result. Note that the lattice mismatch 
between the Ln202 layers and the Fe2As2 layers affects the 
chemical stability of LnFeAsO. 



ing Ln202 layers and negative charges (electron dop- 
ing) in the Fe2As2 layers. Earlier preliminary experi- 
ment showed si gn o f superconductivity below 10 K in 
GdFeAsOi_a;Fa;.[io| Later, Tc was pushed up to 36 K 
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in a polycrystalline s amp le with a nominal composition 
of GdO0.83F0.17FeAs.Elrl Very recently the Tc value was 
increased to 53.5 K in oxygen-deficient GdFeAsO 1 -a; by 
using high-pressure synthesis, [l^ It is of great interest 
whether the Tc can be elevated further in electron-doped 
GdFeAsO systems. 

Up to now, electron doping in the iron-based ox- 
yarsenides was realized through the chemical substitu- 
tion only at oxygen site in Ln202 layers. Because the 
ionic radius of F~ (1.31 A, GN=4) is distinctly smaller 
than that of O^" (1.38 A, GN=4)[li|, F" substitution for 
O^^ in GdFeAsO leads to more severe lattice mismatch 
as mentioned above. In other word, doping F~ (or oxy- 
gen vacancy) in GdFeAsO is particularly difficult, which 
is probably the main obstacle to elevate Tc- Substitution 
of Ln'^+ by relatively large tetravalence ions is an alter- 
native route to introduce electrons. A successful example 
was the electron doping in Ln2-2:Ge2;Cu04 (Ln=Pr, Nd 
or Sm), which has led to the discovery of n-type cuprate 
superconductors^!?]. Th'*+ is a very stable tetravalence 
ions and is as large as Gd^+[lJ|, therefore, we pursued 
the Th^+ substitution for Gd''+ in the GdFeAsO system. 



Figure 2a shows the X-ray diffraction (XRD) patterns 
of the Gdi-j^Thj-FeAsO samples. The XRD peaks of 
the undoped compound can be well indexed based on 
the tetragonal ZrCuSiAs-type structure, indicating single 
phase of GdFeAsO. As for the Th-doped samples, small 
amount of unreacted Th02 can be seen. The refined lat- 
tice parameters are a = 3.9154(2) A and c = 8.4472(4) 
A for the parent compound, basically consistent with the 
previously reported values [l4|. For the Th-doped sample 
of a;=0.2, the fitted lattice parameters are a = 3.9161(2) 
A and c = 8.4386(3) A. The other Th-doped sample of 
x=0.25 has a similar unit cell with a=3. 9166(3) A and 
c = 8.4382(5) A. Therefore, it is evident that the Th- 
doping tends to expand the lattice within a6-planes but 
to shorten the lattice along the c-axis. The substantial 
change in cell constants indicates that thorium is incorpo- 
rated in the lattice. Meanwhile, the a6-plane expansion 
reflects the effective electron doping onto FeAs layers. 

More direct evidence of Th incorporation into the lat- 
tice comes from the chemical composition measurement 
by energy-dispersive X-ray (EDX) microanalysis. Fig- 
ure 2b shows that the microcrystal in the SEM image 
contains remarkably Th in addition to Gd, O, Fe and 
As. Quantitative analysis gives the Gd:Th:Fe:As ratios 
as 0.83:0.19:0.96:1.00 for the a;=0.2 sample (Here we omit 
the oxygen content because the amount of oxygen cannot 
be measured so precisely by EDX technique). This result 
demonstrates that most of the Th was successfully doped 
for the sample of a;=0.2. Obviously, Th-doping level is 
higher than that of F-doping in GdFeAsO lO, illi]. 

Figure 3a shows the p{T) curve for Gdi-^^Th^^FeAsO 
samples. For the undoped GdFeAsO, the p{T) curve 
exhibits an obvious anomaly at 128 K, characterized 
by a resistivity drop with decreasing temperature. In 
LaFeAsO a similar resistivity anomaly was found at 150 
K,[l|, [T^ which has been recently suggested to be asso- 
ciated with a structural phase transition and /or an an- 
tiferromagnetic spin-density-wave transition[i^ [l^ |20, 
[21| . We speculate that the present resistivity anomaly 
in GdFeAsO has a similar physical origin with that 
in LaFeAsO. For x=0.2 and x=0.25, such resistivity 
anomaly disappears, instead, resistivity drops abruptly 
to zero below 55 K, indicating a superconducting tran- 
sition. In addition, the linear temperature-dependence 
of normal-state resistivity near Tc suggests possible non- 
Fermi liquid behavior in the present system, similar to 
that in SmFcAsOi_^F^[i]. 

The magnetic measurement on Gdi_:rTh2:FeAsO sam- 
ples was shown in Figure 3b. In the parent compound, 
the high-temperature magnetic susceptibility (x) follows 
the Curie- Weiss law. The fitted effective magnetic mo- 
ments was 7.96 fiB per formula unit, in good agreement 
with the magnetic moment of free Gd^+ ion. Below 4.2 K, 
X drops sharply, indicating an antiferromagnetic ordering 
of Gd^+ magnetic moments. Note that similar behavior 
was observed in SmFeAsO where the Neel temperature 
was 4.6 K.[22] 

For the Th-doped samples, e.g., x=0.2, the normal- 
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FIG. 2: Phase and compositional analysis of 
Gdi_j;Th2,FeAsO samples. a, Powder XRD of 

Gdi-2,Tha;FeAsO poly crystalline samples. The asterisked 
peaks come from unreacted Th02. b, a representative EDX 
spectrum of a crystalline grain of Gdo.8Tho.2FeAsO sample, 
indicating the incorporation of thorium into the lattice. The 
inset shows the SEM image of the identical sample. The 
marked spot is the position where the EDX spectrum was 
collected. 



state susceptibility roughly obeys Curie- Weiss law. Be- 
low 55 K, X drops steeply to negative values, confirming 
the superconductivity observed above. After subtracting 
the paramagnetic susceptibility of Gd'^+ ions, the volume 
fraction of magnetic shielding at 10 K was estimated to 
be over 50%, indicating bulk superconductivity. The dif- 
ferential xfc curve in the inset of Figure 3b shows that 
the Tc(onset) is over 56 K, consistent with the resistance 
measurement. Therefore, Th-doping in GdFeAsO ele- 
vates the Tc by 20 K {Tc^36 K in GdOo.ssFo.iTFeAsLll]). 
Moreover, the Tc of 56 K is among the highest ever dis- 
covered in iron-based oxypnictides. 

Our observation of superconductivity in 
Gdi_a;Tha;FeAsO indicates that the Ln-site substi- 
tution in LnFeOAs is feasible to realize electron doping, 
hence the high temperature superconductivity. More- 
over, the electron-doping is more easily induced by the 
Th """-substitution compared with the F~ -substitution in 
GdFeOAs. This manifests that the lattice match between 
Ln202 fluorite-type layers and Fe2As2 antifluorite-type 
layers is important not only for the chemical stability of 




50 100 150 200 250 300 

T(K) 



T" =4.2 K 



□-x=0 
•-x=0.2, FC 
A-x=0.2,ZFC 




50 100 150 200 250 300 

T(K) 

FIG. 3: Electrical resistivity (p) and magnetic suscep- 
tibility (x) of Gdi-^ThxFeAsO. a, p versus T. The data 
are normalized to psooA' as the resistivity measured on poly- 
crystalline samples is often higher than the intrinsic value due 
to the grain boundary and surface effect. The inset is an ex- 
panded plot to show the superconducting transition at 56 K 
for the sample of x=0.2. b, x versus T. The lower inset 
shows the differential Xfc curve of Gdo.8Tho.2FeAsO powder 
sample, which locates the onset Tc at 56.3 K. ZFC, zero-field 
cooling; FC, field cooling. 



the parent compounds, but also for the electron-doping 
level in LnFeAsO. It is thus expected that the thorium- 
doping strategy can be applied to other iron-based 
oxypnictides. 
METHODS 

Sample preparation 

Polycrystalline samples of Gdi_a;Tha;FeAsO were syn- 
thesized by solid state reaction in an evacuated quartz 



tube. All the starting materials (Gd, Gd2 03, Th02, Fe 
and As) are with high purity (> 99.95%). First, GdAs 
was presynthesized by reacting Gd tapes with As pieces 
in vacuum at 773 K for 10 hours and then 1173 K for 12 
hours. Similarly, FeAs was prepared by reacting Fe pow- 
ders with As shots at 773 K for 6 hours and then 1030 
K for 12 hours. Then, powders of GdAs, GdaOg, ThOa, 
Fe and FeAs were weighed according to the stoichiomet- 
ric ratio of Gdi_xTha;FeAsO. The weighed powders were 
mixed thoroughly by grinding, and pressed into pellets 
under a pressure of 4000 kg/cm^ in an argon-filled glove 
box. The pressed pellets were wrapped with Ta foils, and 
sealed in an evacuated quartz ampoule. The sealed am- 
poule was slowly heated to 1423 K, holding for 48 hours. 
Finally the samples were rapidly cooled to room temper- 
ature. 

Structural characterizations 

Powder X-ray diffraction was performed at room tem- 
perature using a D/Max-rA diffractometer with Cu Kq, 
radiation and a graphite monochromator. The XRD 
diffractometer system was calibrated using standard Si 
powders. Lattice parameters were refined by a least- 
squares fit using at least 20 XRD peaks. Energy- 
dispersive X-ray (EDX) spectra were obtained by us- 
ing the Phoenix EDAX equipment attached to a field- 
emission scanning electron microscope (SIRION FEI). 
The ground sample powders were placed directly on the 
copper sample holder for making the SEM specimen. 

Electrical and magnetic measurements 

The electrical resistivity was measured with a stan- 
dard four-terminal method. Samples were cut into a thin 
bar with typical size of 4mmx2mmx0.5mm. Gold wires 
were attached onto the samples' newly-abraded surface 
with silver paint. The size of the contact pads leads to 
a total uncertainty in the absolute values of resistivity 
of 10 %. The electrical resistance was measured using a 
steady current of 5 mA, after checking the linear I ~ V 
characteristic. 

Temperature dependence of magnetization was mea- 
sured on a Quantum Design Magnetic Property Mea- 
surement System (MPMS-5). For the measurement of 
the undoped compound, the apphed field was 1000 Oe. 
For the measurement of the Th-doped superconducting 
samples, both the zero-field-cooling and field-cooling pro- 
tocols were employed under the field of 10 Oe. 
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